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1. INTRODUCTION 

Transient synchronous generator models extensively applied in many areas take the part of 
important parts as the beginning point of any power system stability studies. Various generator models have 
been developed [1]—-[4]. For such analysis studies, the models of simple transient states based on the common 
method, such as the primary swing stability of electromechanical oscillation of the rotor, is useful but not 
suitable to envisage generator performance concerning control experiments; it creates dissimilar results with 
the measurements [5]-[8]. 

The hypothesis that the reciprocated inductances among armature, dampers, and field on direct axis 
windings are similar is the basis of conventional stability theory of the generator model. Moreover, both copper 
losses and machine slots are ignored and the latitudinal distribution of the stator fluxes and apertures wave are 
sinusoidal wave forms. The damper winding that are placed near air-gap present flux linking damper circuits 
whose value is nearly equal to flux linking armature. This assumption accelerates satisfactory outcomes for 
many stability studies particularly individuals of network side [9]—[15]. Nevertheless, once it tackles to field 
current explores, there is a substantial error. Concerning transient condition, the measurements of field current 
have exposed greater alternating amplitude [16]. Traditional model of two-axis framework designates only the 
armature circuit appropriately. The additional inductance that indicating the disparity between field-damper and 
field-armature reciprocated inductances on d axis has to be involved in the model [17]. This inductance X¢ 
named Canay inductance, takes significant role to the leakage flux ®-, shown in Figures 1(a) and 1(b) once 
Canay inductance is neglected, together base current of field winding and inductance of armature leakage are 
regarded free parameters. During assessing an identical circuit model, the residual parameters are attained for 
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the purpose of the model replicates same behaviour from terminal viewpoint. The leakage inductance may 
possibly not be equal to the value provided by the real system. One method to creating the rotor circuit is to 
model such effect by a low-order ladder circuit of RL branches in order to look like the complex frequency- 
dependent nature of the rotor’s surface impedance [18]. Another approach is to model it in conjunction with that 
traditionally viewed the partition of current paths from the physical structure [19], [20]. In addition to both 
physical field and damper circuits, a third circuit for the eddy currents induced could also be inserted in the pole 
surface. These three current paths are in near closeness in the rotor slot. 

This study describes the method whereabouts synchronous generator model for stability 
investigation constructed on the hypothesis that the mutual inductances among armature, dampers, and field 
on direct axis windings are not the same. Procedure for making the proposed equivalent model and how the 
advantage of proposed model compare to formerly conventional model are demonstrated. Section 2 tells 
research method. Section 3 describes result and analysis. While, the conclusion is described in last section. 
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Figure 1. Structure of coupling flux components in a rotor slot (a) in theta direction and (b) in radial direction 


2. METHOD 

In order to compute the machine parameters correspondents to the actual measured responses, the 
rotor winding quantities, a better arrangement can be attained with enhancementsts in the structure of the d 
axis model. Shielding impacts of the damper currents and eddy currents induced in rotor can influence 
transient characteristics significantly. The value of this effect depends on the rotor construction. 


2.1. Higher order model of synchronous generator 

The traditional synchronous generator model is obtained from the dealings of the joined stator and 
rotor electric circuit by converting the variables of generator hooked on dg axes. This converting is used 
since the state-space equations depiction is more compact and it enables the analytical method proposed in 
this. The d axis is counted to be leading the g axis by 90 degrees and generator agreement is used to 
described the voltage equations. The traditional generator model is expressed by the corresponding circuits 
shown in Figure 2 with two damper windings both in g axis and in d axis, and counting one field and one 
damper windings, and also stator windings described through q0 and do. 
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Figure 2. The traditional model of g axis and d axis equivalent circuits 


Using relation between magnetic and electric circuits described in [21], [22], each of the flux 
components coupling that the current paths may be transformed to a mutual inductance in an equivalent 
circuit representative of rotor where the identity of these current paths are preserved. Figure 3 shows an 
improved qd circuit model with three g and three d rotor circuits. The inductance, Li.,, and Li.,,, are the 
coupling inductances associated with the flux components in the rotor slot, such as ¢,,, and @;2,. The 
inductances, L},,, Liz, and L}3,, are from the self-leakages of the field, damper, and eddy current. 

Those equation aimed at equivalent circuit through three rotor circuits in Figure 3 [23] can be shown 
as (1). 
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Wq ~ Wg as Xishg 

Wa-Wma = Xisla 

Wo = Xisto 

Whas — Wma = (X3¢ + Xp2cdikaz + Xroctkaz + Xracly 

Wkaz — Wma = Xractkaz + X2e + Xric + Xrecdikaz + Xrac + Xric dig 

Wr — Wma = Xrackeazs + Orie + Xracdikaz + Xte + Xt20 + Xrrc)i¢ 

Wheaqa = Wmq = Xikg3 ikg3 

Wheqa od Wg = Xikg2 ikq2 

Waa iz Wg = Xing tkqt (1) 


The flux linkage equations of the three rotor circuits on the d-axis with mutual coupling can be described 
as (2). 


U r 1, r 1, yf 
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ee 
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Determining parameters of generator model to acceptable the frequency response test data, a correct model is 
decided by first selecting the suitable circuit representation. Consequently, to overlook the leakage couplings 
of three d axis rotor circuits is not analogous to dumping the off-diagonal term of full X matrix in (2), as 
diagonal values from a parameter suitable could be different between without and with off diagonal elements. 

Using the inverse relation B = X77 in [24], ikq3> Ueqa» and if can be determined from the values of 
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The equation for the stator gd currents and rotor q axis currents can be attained in the same way, namely: 


q > Wq [2 Wmg/X1s ig = Wa ~ Wmal/X1s 
ikg3 S (Wigs a Wing) | Xtegs ikg2 = (Wkqo — Wiriia)) Minas ikqi = (Wiegt — Wing) /Xteqa (4) 


the flux linkages, Pq and mq, can be represented in term of the total flux linkages of the winding, that is: 
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The flux linkages of stator and rotor circuit are attained by integrating their individual voltage 
equations: 


Wq =p S{vq — (@,-/@p)Wa + CE Xi) Vina = q)}at 
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Wkaz = —@pTua2 J {bo1Wkar + b22Wka2 + bos + (bar + baz + b23)Wma fat 
Wr = (wy7/Wma) S{E; — Xmd (b31Wka3 + b3iWkaz + b33)p — (bz, + b32 + b33)Wma) fat (6) 


where Ey = Xmap / 1; . The rotor speed, w,.(t), is shaped from expression of slip speed: 
t 
w,(t) We = (P/2J) J, (Ti + Tmech ~ Tate) dt (7) 
where the value of T,,, is negative and it is computed from: 


Tem = (3/2)(P/2w5)(Waig — Wgia) (8) 
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Figure 3. A proposed model with leakage coupling between d-axis circuits [23] 


2.2. Synchronous generator simulation in MATLAB/Simulink 

A higher order of synchronous generator model is designed exposed in Figure 4, by means of 
Simulink components [25]. As shown, the conversion of input stator abc voltages to the rotor gd reference 
frame are made within the abc to qd0 block. The block gd gen comprises simulation of generator suitable in 
rotor’s reference frame. The stator abc voltages, excitation voltage, and applied mechanical torque on the 
rotor are the input of the reference frame. While, the outputs exist terminals of voltage magnitude, current, 
active power and reactive power, and also load angle, slip, and accelerating torque. 

Figures 5 and 6 describe qg axis and d axis circuit parts within the gd gen block that contains three 
rotor circuits in every axis. The block of g axis circuit comprises three damper windings. The d axis circuit 
encloses one field armature winding and two damper windings. The synchronous generator parameters are 
shown in Table 1. 


Table 1. High order generator parameters [26] 


Machine type-1 Machine type-2 
Characteristic : ‘ Characteristic : 5 
Constant q axis part d axis part Constant q axis part d axis part 
Frated 60 xmd 1.70 xmq 1.61 Frated 60 xmd 1.70 xmq 1.61 
Poles 2 xpr2c 0.06523 xplkq3 0.1225 Poles 2, xpr2c 0.06650 = xplkq3 0.1225 


Vrated 26e3  xpric 0.01925 xplkq2 0.3248 Vrated 26e3  xpric 0.0 xplkq2 0.3248 
Srated 722e6 xple 0.1055 xplkq1 0.6802 Srated 722e6 xple 0.12094 = xplkq1 0.6802 


rs 0.004 = xp2c 0.03076 — rpkq3 0.2237 rs 0.004 = xp2c 0.06563 = rpkq3 0.2237 
xls 0.19 xp3c -0.009134 = rpkq2 0.03537 xis 0.19  xp3c 0.00988 = rpkq2 0.03537 
xpd 0.346 = rpkd3 0.02467 ~— rpkq1 0.005698 xpd 0.346 rpkd3 0.02334 _—s rpkq1 0.005698 
xpq 0.642 = rpkd2 0.01297 xpq 0.642 rpkd2 0.01848 

H s) Ipkd! 0.00112 H 3 Ipkd1 0.00109 


Higher order model of synchronous generator (Sugiarto Kadiman) 


1446 O 


ISSN: 2088-8694 


Workspace 
Mux Mux 
F e t, 7 7 
omega’t vl i 987 gen Power ‘tem —Lipf wriwb | ia ib ic ikd ——Jikg 
arfgle 
Mux Mux Mux d 
? 2 +) Clock VIPQ qdr2abc 
e 
Mux3 Vm3 Mux2 V2 Mux! Vmi 
a . 
—— iq_gen = 
>| _ulereostutth) vag 4 | “ah Ground 
Fen vbg i 
SS veg vd wriwb 4 
{__] ypyrcosutt}-T2pibys) Tem 
Font we inf mani 
ie id_gen osc 
>| ulzrrcos(ult}+T2piby3) eT 4 
Fon2 
gad _gen sin (wrt) 
£7 cos (wrt) 
Ef 
Tmech = 
Figure 4. Simulink diagram of high order synchronous generator model 
4s ‘s 4 Out_psiq 
wb*(u[2]-u[3]+(rs/xls)*(u[1]-u[4])) s + 
Fen psiq L-+| f(u) T >( 3 ) 
In_wrpsid Fon3 Out_psimq 
Mux 
Mux L-+| wb*rpkq3*(u[1]-u[2])/xplkq3 tL -—>| -(u[1]}-u[2])/xIs |/——______+(2) 
Fen2 psipkq3 Foen4 Out_-iq 
Mux2 
um 
+ 
Mux >| a y 1 
wb*rpkq2*(u[1]-u[2])/xplkq2 si 
Fen5 psipkq2 
Mux5. 
Mux t-+| wb*rpkq1*(ul1]-ul2])/xplkq1 i 4 
Fon6 psipkq1 
Mux6 


Or 
@— 


In_wrpsiq 


~) 


wb*(u[2]+u[3]+(rs/xls)*(u[1]-u[4])) 


Fen 


Out_psid 


(u[1]-u[2])/xIs. 


xMD*(u[1/xls+u[2]*b1col+u[3]*b2col+u[4]*b3col) HG) 
Mux Fen3 Ducesind 
| Mux3 
[ 
Loe -wb*rpkd3*u[1] t | 
Fen2 psipkd3 
L_| -wb*rpkd2*u[1] t + 
Fon6 psipkd2 
a Demux Out_ipf 
wb*rpf*(u[1]/xmd-u[2]) >| i + I. | Muxd 
Fi Sum2 
Font psipf 
Mux1 


Figure 6. The inside part of d axis block 
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3. RESULTS AND DISCUSSION 

With the resolve of exploring shielding effects of the damper winding currents and the eddy currents 
produced in rotor can impact transient characteristic of the rotor, numeral simulations were conducted for 
various different operating conditions, for example step change in mechanical torque and step change in input 
voltages. The high order model of synchronous generator is used and the values of field, torque and input 
voltage perturbations about operating point of field current are subject to change in p.u. Figures 7(a)-7(c) 
show curves of field current for the period of short circuit simulation is at filed, torque and input voltage 
perturbations using machine parameters in Table 1. Field current under a three-phase sudden short-circuit is 
described by (9) [27]. 


y 
ip = ipo + ifo (4) fet/t4 - (1 - wie) en t/Ta — moe cos anf} (9) 
Where if, is field current before abrupt short circuit; while Tq denotes time constant of d axis damper 
circuit, and others are general reactance and time constants expressed by the two-axis theory. The magnitude 
of the field current will vary with the existence of perturbations. The emergence of torsional disturbances 
causes the value of the field current to decrease. Conversely, a disturbance in the input voltage will increase 
the field current. but after the disturbance subsides the field current value returns to its previous value. 
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Figure 7. Field current curve under three phase short circuit on: (a) field perturbation, (b) torque perturbation, 
and (c) input voltage perturbation 


Figure 8(a) shows curve of stator input voltage. The form of the signal is a step function with a 
magnitude of | p.u; in the period between 0.1-2.7 seconds has a magnitude of 0 p.u which is considered as a 
disturbance. Figures 8(b) and 8(c) show curves of field current during short circuit simulation under different 
two of machine parameters. The use of a larger reactance value in X,;, and X,2, results in oscillations that 
occur in the event of a disturbance. 
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Figure 8. Short circuit simulation of a proposed model with leakage coupling between d-axis circuits 
(a) stator input voltage, (b) short circuit response from 2X3 model with non-zero x,,, and x,2, with machine 
type-1, and (c) short circuit response from 23 model with non-zero x,1, and x;2- with machine type-2 


4. CONCLUSION 

The authors examined synchronous generator dynamics, for example field current characteristics at 
transients, three-phase sudden short-circuit in specific. The following goals are described: i) To observe the 
field current effect based on the quantity in the three-phase abrupt short-circuit experiment, a model of 
synchronous generator concerning field mutual leakage reactance, namely Canay reactance, that relates to 
magnetic flux connecting merely on the rotor area should desire be applied; ii) The simulation results display 
that the field current in the short circuit assessment will be dissimilar due to occurrence of different types of 
perturbations; and iii) Considering dynamic stability, the resulting trends are noticed in filed current 
fluctuating frequency for application of high excitation: a. designed for x,, > 0, negative damping trend 
rises; b. designed for x,, < 0, negative damping trend recovers. 
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